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ABSTRACT 


The photodeeorbtlon of molecules and atoms from the surfaces of inter- 
stellar grains can be an Iji^ortant source of heating for the interstellar 
medium and the origin of instabilities which may separate grains and gas. 

For low densities, the force exerted on the grains is proportional to the 
gas density and independent of the radiation intensity; for high densities, 
it is proportional to the radiative flux and independent of the gas density. 
This force may act differently on grains of different sizes. The photo- 
electric effect may also be an efficient mechanism for the separation of 
gas and dust in diffuse clouds. 
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I, Introduction . 

The non-lonizinc radiation from stars has an energy density comparable to 
that of the other con^onents of the interstellar medium, but its effects on the 
overall structure of HI regions have generally been thought to be small. We 
point out that this may not be *the case, and that starlight, through its inter- 
action with atoms or molecxiles adsorbed on the surfaces of grains, may play an 
important role in many parts of the interstellar medium. In this note, ve 
present simple physical arguments which suggest the circumstances under which 
the effects of starlight may be important, and which will enable us to estimate 
the contribution of this process to the heating and the stability of the inter- 
stellar medium. 

In Section II we discuss how the photodesorption of atoms or molecules 
adsorbed on the surface of interstellar grains may be a mechanism for the 
strong coupling of the grains to the stellar radiation field. In Section III 
we investigate the role which, this process plays in the separation of gas and 
dust in the interstellar medium, and show that the canonical assumption of a 
constant ratio of gas to dust, while it may be true in some average sense, 
cannot be true for individual clouds. We shall exhibit an instability which 

f 

occurs in uniform gas-dust clouds which spontaneously separates the tviro, even 
In the presence of a uniform stellar background. In Section IV v/e sliow that 
this process may also be an important heating mechanism for diffuse interstellar 
clouds which is not affected by the attenuation of radiation in the cloud. 
Finally, in Section V the different consequences of these results arc 
briefly discussed. 
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II. Photodesorption and the Force on Single Grains . 

Our discussion begins with the recognition by Reddish (1971) that, if 
near UV photons can liberate adsorbed molecules from the surfaces of inter- 
stellar grains, the momentum transferred to the grains will be enhanced by 
a large factor over that due to radiation pressure alone. The reader may 
observe a similar effect in a device usually known as the Crookes radiometer, 
which revolves in the opposite sense fron that which one would expect frcan 
radiation pressure. (A Crookes radiometer in an extremely high vacuum goes 
around the "right" way. ) 

Although the photodesorption of atoms or molecules from solid surfaces 
has been Imown experimentally for a long time, it still is not clear which 
are the basic mechanisms responsible for it. Watson and Salpeter (1972) 
argue that the adsorption of a photon involves the transition to an excited 
electronic state of the adsorbed atom or molecule. By the Franck-Condon 
principle the transition occurs at a constant value of the intemuclear 
distance (i.e. the distance from the surface to the adsorbed particle). 

Since the adsorption potential for the excited state is likely to be quite 
different from that of the ground state, the difference will usually be 
greater than the adsorption energy. If the particle proceeds to a distance 
of the order of several Angstroms before returning to the electronic ground 
state, it will leave the surface. 

For sites of enhanced binding, on which particles arc chemisorbed, 
however, the electron may be excited into a state which is equivalent to 
one of the electronic states in the solid. If the overlap between this state 
and the states inside the solid is large, the probability of desorption e 
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wlll be small; if the excited state Is localized, c may be large even for 
chemisorbed species. For highly dissimilar species boxmd to the crystal, 
the overlap Is smaller, because of the difference in the adsorption bond 
length and the lattice constant. 

The argument given above also leads us to assume that most of the photon 
energy is given to the desorbed particle, since the energy transfer to the 
lattice Is expected to be small. If the desorbed particle leaves the surface 
in. an excited electronic state, this state will be of a well-defined energy, 
and the difference between the photon energy and that of the desorbed excited 
state will also be precisely defined, so the difference between the photon 
energy and that of the excited state will appear in the only available form, 
as kinetic energy of the particle. We are therefore led to assume that a 
large fraction of the photon energy appears as kinetic energy of the desorbed 
particle. Other mechanisms that have been proposed for photodesorption are 
listed by Greenberg (1973). Recently, Gersten, Janow and Tzoar (1975) 
have made a detailed theoretical study of the photodesorption of single 
atoms. They estimate a quantum yield for photodesorption from a mono- 
layer on a surface of ~ 2 x 10 in agreement with some earlier experimental 
work by Fabel, Cox and Lichtman (1973). Of particular interest here is the 
possibility, mentioned by Greenberg, that the molecule goes up to on excited 
electronic state, returns to the ground electronic state, but in an excited 
vibration-rotation level. This internal kinetic energy may be converted to 
translational energy away from the surface by collision with surface atoms. 

In this case the energy of the ejected molecule is related to the vibrational 
levels of the molecule rather than to the energy of the photon. Present 
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experiments do not allow to discriminate among the different possibilities; 
Oreenberg, however, from experimental data on several molecular species was 
able to estimate that the relative quantum efficiency is very near one, in 
accord with the arguments of Watson and Salpeter, so that an appreciable 
fraction of photons absorbed -lead to desorption of molecules . The force 
exerted on a grain depends on the square root of the kinetic energy of the 
molecule; therefore, at this stage, we assume that the kinetic #»nergy 
of the ejected molecule is essentially that of the photon. 

A grain in which a fraction f of the surface area a is covered by 
adsorbed atoms or molecules of mass m, experience a force 

F " Bo f m^ e \ H (hv)^ dv , (l) 

where is the flux of starlight energy and c the probability that n photon 
will release a particle from a fully covered surface. This estimate is under 
the assumption that the particles desorbed leave the grain isotropically; it 
only changes slightly for other angular distributions, B is an efficiency 
factor, near unity, which includes the effect of diffraction. If we take 
f *s 1, e = 1, and m to be the mass of a hydrogen molecule, this force is 
about 3 X. 10^* times stronger than the radiation pressure. It is therefore 
apparent that this grain Jet effect, as we may call it, may be important if 
e is large even if only a smaH fraction of the grain is covered by adsorbed 
atoms or molecules, or, vice versa, with a smaller e on a completely covered grain. 

The photodesorbed particle and the grain carry equal and opposite 
momentum; the net momentum transferred to the gas by tliis process is tlms 
aero, if we neglect the momentum of the initial photon. The grain Jet effect 


thus causes motion of the grains through the gas; the radiation exerts a force 
much larger than the radiation pressure on a small fraction of the mass of the 
gas-dust system. This means that the instabilities we shall find as a result 
of this force will grow very fast compared to the dynamical timescale associated 
with the gas itself. 

In order to estimate f, we must consider the processes by which adsorbed 
particles stick to and leave the surface of the grains. An extensive discussion 
of these processes has been given by Hollenbach and Salpeter (1971, hereafter 
referred to as HS). They cons! ^ered the process by which molecules stick to 
the grain, and the rate of thermal evaporation, but not photodesorption. 

The average time for a particle to evaporate from the surface of a grain 

is t * (l/v ) exp(+D/T„) in the notation of HS, where D is the adsorption 
ev o g 

energy, T the grain temperature and v the characteristic vibration frequency 
g 0 

of the adsorbed particle. The grain temperatures and the adsorption energies, 
of course, vary greatly according to the nature of the grains. For example, 
the adsorption energy of H atoms on a perfect ice surface was estimated by HS 
to be about 450® K, and for H^ about 550® K. They estimated similar adsorption 
energies for graphite surfaces. However, binding energies of both atoms and 
molecules to defect sites on the grain siirfaces can be manj' times this energy, 
up to chemisorption energies of several electron volts on highly reactive sites. 
If particles are adsorbed, or chemisorbed, on such enhanced sites, they will 
remain so for all reasonable grain temperatures, thus setting a lower limit to 
the value of f. Unfortunately, we do not know how to estimate the number of 
enlianced sites vdiich may bo present on a groin; we therefore shall nol incliidc 
the contribution of ouch sites in our estimates of f. It must be rememherod 
that, if tlicre is on appreciable number of such enhanced sites on a typical 
grain, the occupution number, f, will be larger and the groin jot effect 
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correspondingly stronger. 

The rate at which molecules stick to the surface of the grain Is 

t ■ (o V n S), idiere n Is the gas density and S is a sticking probability, 
0 

estimated by HS to be about 0.9. v is a sooeirtiat complicated function of 

the gas-dust velocity, but has the simple limiting foxms v «■ c , the Iso- 

0 

thermal sound speed, If the dust Is nearly stationary with respect to the 
gas, and v ■ i^gas-dust supersonic grain velocities. Th , rate at which 
grains evaporate Is f t'^ and the rate at which they are photodesorbed is 


A* J 

t^^ * f e • Aira J dv s f e o ^ . 


( 2 ) 


We may determine f by equating the rate of sticking to the 6\sn of the 

evaporation and photodesorption rates. As was shown by HS, the evaporation 

rate becomes large if T exceeds 

g 


T •» D/ln(t V ) 15° K . (3) 

c so 

In that case, f will be negligibly small, and we may neglect the process under 
discussion. Estimates of free-space grain temperatures (e.g. Leung, 1975) give 
temperatures for silicate grains of 13° K, for graphite grains of 30° K, for 
silicate-ice grains of 13° K and for silicone carbide grains of 48° K. It is 
thus apparent that, for at least some types of grains which have been proposed, 
we may neglect evaporation under typical interstellar conditions. 

If we neglect evaporation from the grains, then we may estimate f as singly 


f ■ (S/e) n v/^ . 


(0 
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Thl8 forn Is, of course, applicable only If f « 1. If f is greater than one 
by this prescription, we shall set it equal to one for purposes of our estimates, 
f is nearly one for densities greater than n^, the critical density. 


■ (♦/v) it/s) . (5) 

This density separates the regimes in idiich the grain Jet effect has qualitatively 
different depend^ce on gas density and on the radiation field. 

If n < n^, the force on eaa;' grain is given 

■ n V o S< (hv)^> • 6 , (6) 

where 6, the anisotropy of the radiation field, is defined as 


« - 


_l 


t H/(h'))*dv 

So 


(7) 


This force has several remarkable properties. It is independent of the intensity 
of the radiation field, and depends only on its degree of anisotropy. This be- 
havior is due to the assus^tion that evaporation can be neglected; in that case, 
f is inversely proportional to the flux of photons, so the number of photode- 
sorptions per second is constant (and equal to the number of atoms sticking to 
the grain). For the same reason, the force is independent of c, and proportional 
to n. 

The drag force on a grain (moving through neutral goo) is also proportional 
to n, and to the square of the gas-grain velocity for supersonic motion. Thus, 
the velocity of the groins with respect to the gas is approximately 1/2tt times 
the velocity with which the particles are ejected frean the grain; for typical 
interstellar photon energies of 6 c.v. this corresponds to v ■ 8 km/sec. • 6. 
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in the opposite Halting case, n » n^, the force la 

F ■ o 6 c f J (hv)^ dv (8) 

.proportional to the Intensity of the radiation and to e, and Independent of 

n. This is qualitatively similar to the ordinary radiation pressure, but 

larger by a factor c(idc /hv), fdiere v is the average photon frequency. 

In the solar neighborhood the average energy density of radiation Is 0.01 

photons cn*'^, Witt and Johnson (1973), and if ve assume a mean photon energy of 

5 -1 

6 e.v. and a velocity v of about 10 cm sec , the critical density n has a 

c 

value of lO^ ^cm*'^ e . According to the estimate of Gersten et al (1975), 
e > 10”^ for single atoms. It appears, therefore, that for tlwse diffuse clouds 

m0 

in the solar neighborhood where the grain temperature is low enough, T < 13 X, 
the radiation exerts a force on the grains independent of the radiation Intensity 
and proportional to the gas density. Denser clouds are in the second regime. 

It is important to remember that the grain Jet force at optical depth 
T • i ln(mc /hv) 10 + In c can be comparable to ordinary radiation pressure 
at T ■ 0. 

III. Separation of Gas and Dust. 

The grain Jet effect is the source of Instabilities in the interstellar 
medium which may lead to the separation of gas and dust and to the separation of 
grains of differing composition and/or sizes. If we consider a uniform distribu- 
Xion of gas and dust immersed in a homogeneous and isotropic radiation field, 
and subjected to a perturbation in duct density alone, such o perturbation is 
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unatable. . Hie origin of this inatobility la phyaloally obvious; a' region of 

inereaaed density attenuates the interstellar radiation field, inroduclng an 

* 

anisotropy in the intensity. This anisotropy in turn drives the grains toward 
the region of higher density, amplifying the perturbation. 

It is a short calculation to demonstrate that the system described above 
is unstable to perturbations in the density of the dust component. If we con- 
sider the gas and dust to be two fluids interacting only through the drag force 
on the moving grains, we may XK>te that the grain jet effect does not transfer any 
net BKSDentum to the gas, except for that given by the radiatim pressure on the 
grains, i^ch we shall neglect. In order to be consistent, we therefore neglect 
all terms of order 1/c in the equatic^ of motion of the dust. This is equivalent 
to neglecting the difference in the radiation field as seen in the frame moving 
with the dust from that seen in the frame fixed to the stars. We may then write 
the equations of motion for the dust in the form 


♦ V(pv) ■ 0 

-W ■ - B ”h * A V'O 

where the coefficients A and B are given by 



(9) 

( 10 ) 


( 11 ) 


f« 1 


The only intone of interest in our prohlen are the ultraviolet ones 


oapable of desorbing '^toeis and/or nolecules froa the grains; we any thus 
neglect *Uie eoiisslon of the grains and write the transfer equation as 

y^-.kl + ko J* , (13) 

where a is the albedo of the grains <uid is the ^nisslon by the \mifonn 
background of stars. These eq\»tions are the starting point for our linear 
stability analysis. 

We introduce a perturbation in density, * p exp [nt - ikZ], rlong 
with corresponding perturbations in velocity and radiative intensity. The 
linearized equations of notion and the equation of transfer become 

npj^ ■ 1 k Vj^ pQ (14) 

nVj^ * - B n^ V ♦ A H^/Jq (15) 


Here, • k p^ and y ** cos 6 as usual. The equation of transfer applies to 
one frequency, but, since it is linear, we may regard it as applying to all 
frequencies of interest, with the appropriate opacity. 

To sinplify the solution of this system for n(k), wo introduce the 
dimensionless variables ic - k/iCQ,* n(n - B Uj^); 

eliminate p^^ and v^^ by using (14) and (15), 
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it obUla 


(1 - 1 K \>yt^ ♦ 1(1 - a) Cl ■ “ Hi 


(17) 


Integrtiing (17) over v, «• have 


»1 - i K ® , 

while, If we Mlve (17) for ax^ then integrate, we obtain 


(18) 


- -jf [a III - 1(1 - «) A 

♦ 


(19) 


ahera 


1-1 tan^^ic 


Iy T 


i t 

~1lW 


( 20 ) 


Conbining the relations (18) and (19)> we ni^ ellninate and finally obtain 
the result 


7 “ O - it 


( 21 ) 


In order to put this result back into dimensional form, we introduce 
the characteristic time * (k^A)“^ , so that 


"<" * V ■ "o ^ • 


( 22 ) 
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For optically thin perturbations (k » 1), we may neglect the second term and 
slowly obtain 

n(n B njj) ■ (l - a) . (23) 

Notice that. If the albedo of the grains is 1, no photons are absorbed and 

2 

there is no photodesorption; it is therefore satisfying .to note that ^ ■ 0 
in this case. 

Before discussing the application of this result, we may give a short 
heiiristic derivation of the growth time Tq. If we consider a perturbation of 
wavelength X (■ 2“ir/k), the change in optical depth in a wavelength is about 
^1 * The 2 *efore, the anisotropy of the radiation field will be 

approximately d * grains, wliich quickly achieve their limiting 

^ed of A/B are driven through the gas by this anisotropic radiation. 

In order to increase the magnitude of the perturbation by an amoxint comparable 
to 6p, the grains must move a distance L * X/6. At their limiting velocity, 

B Ujj » n in (22), 

corresponding to the fact that the grains achieve their limiting velocity 
very rapidly, we then have n ^ (B n^jA^A)*^! - a + (1 - cx)^/(ct - iy)} . 

Since the factor in the cxirly brackets is just a number near one, the growth 
time derived from this argument is the same as that obtained from the 
stability analysis. 

For the low density regime, which includes most diffuse interstellar 

-9 2 — 

clouds, we may estimate by taking a'V'lO cm; hv»*3e'V.; S~l; 

10 7 

Cg * 1 Km/sec. and “ 10 nijj. In this case, ■ 10 years/nj^. This 

time is less than the typical lifetime of interstellar clouds; we may therefore 

expect the grain jet effect to be of importance in the structure of such clouds. 


this will take a time T = L/v = B If we put 
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The differences in temperature among grains of dlffer^t sizes and 
coDQMsltlon, leading to different thermal evaporation rates (c.f. Eq. (3) ), 
generator a separation among grains with different properties through the 
grain Jet effect. The small grains wlilch radiate inefficiently In the 
Infrared, may become too hot to be driven through the gas by this 
mechanism, while the larger ones may move rapidly with respect to them 
and the gas. The qiiantitatlve details of this process will depend, of 
course, on the (unloiown} composition of the grains themselves. 

In the regions near hot stars, the radiation field Itself may be 
higlily anisotropic, and, as Reddish originally suggested (see also Brand, 

1975), may drive the dust at high velocities. However, the stellar 
radiation may heat the grains enough to make the thermal evaporation 
proceed more rapidly than photodesorption, since the former depends ex- 
ponentially on the grain ten5)erat\ire. The growth rate n of perturbations 
in the low density regime of the grain jet effect does not depend on the 
redlatlon flxix; when the density is high enough, the second limiting case 
obtains and because of its dependence on the flux of radiation, the in- 
stability eventually shuts off. At this stage, hov/ever, the dust density 
should have Increased by a very large factor, and our linearized analysis 
no longer applies. 

The second regime of the force is analogous to the radiation pressure. 

If e B 10**^, the force due. to photodesorption is equal to that due to radiation 
pressure. For those regions in the second regime of the grain Jet effect, 
using the same values of the parameters as before we find » (5 x 10 /c)o years. 


i 
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vhere t Is the opacity. If e 5 lO"^, the growth time again is less than the 
typical lifetime of interstellar diffuse clouds. 

The second regime also applies to the dynamical effects of the photo- 
ejection of electrons from the surface (Watson, 1972). The momentum transferred 
to the grain per photoelectron is smaller than the one transferred by photo- 
desorbed molecules by a factor (m^/m)^, where m is the mass of a typical 
molecule. This effect is con5>ensated, however, by the very high efficiency 
expected for the photoelectric effect. Values of this efficiency are very 
uncertain at present time but according to Watson they range from 0.01 to 0.1. 

For very small grains. Jura (1975) estimates that the efficiency could be as large 
as 0.67. For optically thin clouds, therefore, * 3 x 10^ ~ 3 x 10^ years. 

The point to be emphasized then is that the complete mixtxire of gas and 
dust in an isotropic radiation field is an unstable ccmfiguration: all diffuse 
interstellar clouds should have an inhomogeneous dust distribution. 

According to the treatment of Spitzer (1968) of the equilibrium electric 
charge of grains, the neutral subsonic drag force, used in the preceding 
calculations, dominates over the plasma drag if the ionization degree of the 
gas is. ^ 3 3C lO”^. For HI clouds we are then Justified in neglecting the 
plasma drag. 

In the intercloud medium, hi.'wever, the reverse is true and in this case 
we have to multiply (12) by a factor ~ 363. x f where x 'the degree of ioni- 
zation. Following the conventional view of the Intercloud medium, e.g. Dalgarno 
and McCray (1972), we adopt x ~ 0.1 and n^ " 0.3. In this case ~ 10^ 
years, a value much larger than many other characteristic times of tlie intercloud 
medium. Neither the groin Jet effect nor the photoelectric effect are able to 
separate the dust from the gas in the Intercloud medium. 
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IV. Heating of Interstellar Clouds. 

At the present time the basic mechanisms respcmslble for the heating of the 
Interstellar diffuse clouds are not fully understood. In the low density regime 
of the grain Jet effect » which includes diffuse clouds, the rate of ejection of 
molec\iles or atons Is Independent of the Intensity of radiation and the efficiency 
for photodesorption, being equal simply to the rate of arrival of atoms to the 
surface.. As the adsorption energies are very small compared to the average energy 
of the photons, every atom that sticks to the surface of the grains Is photo- 
desorbed and carries off an energy of several electron volts. This process may, 
tiierefore, constitute an Isqportant heat source for the Interstellar medium, c.f. 
Silk (1973). 

-12 

If we assume an average dust to gas ratio of 10 gralns/H atom, the 
heating rate Is approximately 

H r 5 . njj^(T/80)^ ergs s"^. (24) 

According to Glassgold and Langer (1974), two of the most In^rtant 
heating mechanisms for clouds are cosmic rays and the photoelectric effect 
from grains. The ratio of the heating rate due to photodesorption, as given 
(24)> and the heating rate due to the photoelectric effect, is 
(njj/25) (T/80)^ (y/0.1)~^, where y is the efficiency of the photoelectric 
effect. It appears then, that for diffuse clouds, photodesorption may be more 
important than or at least comparable to the photoelectric effect as a source of 
heating. 

Until Ujj the heating is nearly independent of the radiative intensity, 

and thus Is unaffected by shielding due to the grains. The fact that the heating 
rate has the same density dependence as the cooling rate ought to enhance the 
thermal stability of the clouds. 
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Although there Is no observational evidence for H 2 molecules with kinetic 

energies as high as 6 e*v. In the diffuse clouds observed by the Copernicus 

satellite (Spltzer and Jenkins, 197$), the thennallzatlon time for such molecules 

.5 

Is sufficiently short If n > 10 cm that this Is cosQ>atlble with grain Jet 
heating for such clouds. If the velocities of ejection are smaller It Is 
doubtful that this process Is an In^rtant source of heating; Its dynamical 
effects, however; can still be substantial and this Is the point we want to 
tfBphaslze In this note. Equating the heating rate given by (24) to the cooling 
rate due to CII, the main cooling mechanism In diffuse clouds (Dalgamo and 
McCray, 1972), we estimate a cloud ten 5 >erature of 'V/ 90° K, Independent of 
the density. 

Y. Discussion. 

On the basis of the preceding results it is ten 5 >ting to speculate about 
typical consequences. The cltimps of dust formed by the instability may be 
sites of very rapid formation of molecular hydrogen if they become optically 
thick. Although it is not possible to calculate the local enhancement of the 
dust to gas ratio from the linear theory of instability, simple numerical 
estimates indicate that increases on this ratio of very large factors are to 
be expected. Diffuse clouds should thus have an inhomogeneous molecular 
distribution. 

Because of the heating mechanism, the gas temperature inside these clumps 
of high dust density would increase; observations of molecular hydrogen, such as 
those by the Copernicus satellite, would tend then to sample the regions of highest 
teraperatxire in a cloud. The temperature inhomogeneities would also generate 
motions of the gas to equalize the pressure differences. 

In summary, we may conclude that the grain Jet effect may have important 
consequences in all the neutral parts of the interstellar medium except the 
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dlffuse Inter cloud regions. In the moderate density regions (n < n^)i the 
force per grain Is proportional to the gas density and the radiation anisotropy, 
and Independent of the radiative intensity. The contribution of photodesorption 
to the energy balance of the Interstellar gas may be Important In this regime. 

In regions of high density, the force on the grains Is Independent of the 
density and proportional to the radiative flux. In such regions. It exceeds 
the radiation pressure by large factors. 

If the parameters discussed above have values udilch lie well within the 
range of those usually assumed, this process may be an inqportant source of 
l»at in dense HI regions, the origin of new instabilities whose consequence 
for the formation of dense clouds, molecules and stars must be explored, and 
a mechanism for producing large variations in the ratio of gas to dust In the 
galaxy. 

As with all mechanisms which depend on properties of the Interstellar 
grains, many uncertainties still remain In the basic physics. We hope that the 
discussions in this paper will be a stiimilus for further laboratory studies on 
photodesorption from solid surfaces as well as on research on its astrophysical 
consequences. 

We are greatly indebted to a referee for his constructive comments on 
the original version end to Profs. A. E. Glassgold and F. Kahn, for interesting 
discussions. This research has been supported in part by grant NGR-33-016-196 
from the National Aeronautics and Space Administration. 
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